& Not Peer Reviewed 


bio- protocol Preprint @ Ask a question online 


Cell cluster annotation, trajectory inference, and signature 
enrichment 


e Ekaterini Maria Lyras a) Alexander Mildner 
Updated date: Apr 5, 2023 


j=] An abbreviated version of this protocol was published in eLIFE 
Tongue immune compartment analysis reveals spatial macrophage heterogeneity 


DOI: 10.7554/eLife.77490 


Detailed protocol 


GSVA was performed according to the published protocol (https://bioconductor.org/packages/devel/bioc/vignettes/GSV A/inst/doc/GSVA.html). 


library("Seurat") 
library("GSVA") 
# load your seurat object and the signature genes you are interested in, here: "your.seurat.object" and "signature.genes.1" 
signature.genes.1<-c("Adam19", "Amica1","Ap1s3","Ass1","Bcl11a", "Btla", "Ccr7", "Fit3", "Gpr114", 
"Gpr132", "Gpr68", "Gpr82","H2-Eb2", "Hmgn3", "Kit", "Klri1","Kmo", "P2ry10", "Pvrl1", 
"Rab30", "Septin6", "Slamf7", "Traf1", "Zbtb46") 


gs <- list(sig1=signature.genes.1) 
Y<-GetAssayData(object=your.seurat.object, assay = "RNA", slot="data") 
gsva.es <- gsva(Y, gs, verbose=FALSE) 
your.seurat.object[["sig1"]]<-gsva.es["sig1",] 
FeaturePlot(your.seurat.object, "sigi", cols=c("blue", "lightgrey", "red") ) 
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